Stable two-channel Kondo fixed point of an SU(3) quantum defect in a metal: 
renormalization group analysis and conductance spikes 
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We propose a physical realization of the two-channel Kondo (2CK) effect, where a dynamical 
defect in a metal has a unique ground state and twofold degenerate excited states. In a wide range 
of parameters the interactions with the electrons renormalize the excited doublet downward below 
the bare defect ground state, thus stabilizing the 2CK fixed point. In addition to the Kondo tem- 
perature Tk this system exhibits another low-energy scale, associated with ground-to-excited-state 
transitions, which can be exponentially smaller than Tk- Using the perturbative nonequilibrium 
renormalization group we demonstrate that this can provide the long-sought explanation of the 
sharp conductance spikes observed by Ralph and Buhrman in ultrasmall metallic point contacts. 
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The two-channel Kondo (2CK) effect has been intrigu- 
ing physicists ever since it was proposed by Nozieres and 
Blandin in 1980 [ij because of its exotic non-Fermi liq- 
uid ground state with a nonvanishing zero-point entropy 
It arises when a discrete, degen- 
erate quantum degree of freedom is exchange-coupled 
to two conserved conduction electron continua or chan- 
nels in a symmetrical way Q. Early on, Vladar and 
Zawadowski Q put forward two-level systems (TLS), 
e.g. an atom in a double well potential embedded in 
a metal, as a physical realization of a 2CK system. 
However, this model encountered difficulties because for 
small exchange coupling, J <C 1/A^(0), with A^(0) the 
conduction density of states at the Fermi energy, the 
renormalized level splitting of the TLS cannot become 
less than the exponentially small Kondo temperature, 
Tk ~ exp[— 1/(4A^(0) J)], and a large effective coupling, 
N{0)J « 1, cannot be realized even by invoking enhanced 
electron assisted tunneling via virtually excited states of 
the TLS because of cancellation 6] and screening [7] 
effects. It remains difficult to reach by a resonance en- 
hancement of iV(0) proposed in Ref. [8|. 

Signatures of the 2CK effect have been observed ex- 
perimentally in heavy fermion compounds Q and most 
recently in especially designed quantum dot systems [lB| . 
However, perhaps the most intense interest has been pro- 
voked by the conductance anomalies observed by Ralph 
and Buhrman [llj in ultrasmall metallic point contacts. 
This is because, on the one hand, these measurements 
showed most clear-cut 2CK signatures near zero bias, in- 
cluding the theoretically expected 12, scaling behav- 
ior |14l |. and, on the other hand, these zero-bias anomalies 
(ZBA) occurred in seemingly simple Cu point contacts. 
In addition, sharp spikes were observed in the conduc- 
tance at elevated bias [lH, such that contacts with con- 
ductance spikes always exhibited also the ZBA. These 
results have remained essentially not understood to the 



a) ^ Ji 

W a a // rn=1 t 

XV-IO oy/ J 


^' <:;;o 




<:jl> 


m=0 ' 





FIG. 1: a) Level scheme of the three-state defect, defining the 
coupling constants in Eq. ([ij; b) 2nd and c) 3rd order con- 
tributions to the current through a point contact with Kondo 
impurities. — : conduction electron, : defect state propa- 
gators. A • indicates a sum over all 4-point vertices of Eq. fl}. 



present day mainly due to the lack of a microscopic model 
with a stable 2CK fixed point, although some of the fea- 
tures could alternatively be accounted for by a distribu- 
tion of TLS without invoking the 2CK effect 



15|. How- 



ever, a reliable understanding requires an explanation of 
all the experimental features. The difficulty resides in 
the fact that in the usual 2CK models the Kondo SU(2) 
symmetry is easily destroyed by relevant perturbations. 

In this Letter we present a realistic, microscopic model 
for the 2CK effect, where the Kondo degree of freedom 
are the parity-degenerate rotational states (m = ±1) of 
an atomic three-level system (3LS) and the channel de- 
gree of freedom is the magnetic conduction electron spin. 
Such a level scheme (see Fig. [T|) is generically realized by 
an atom (hydrogen ion) moving in a modulated Mexi- 
can hat potential which, e.g., is formed in the interstitial 
space of the (111) plane of a Cu lattice. Although similar 
models have been considered before [H, [l7[, they have 
not been analyzed in detail, especially not in a current- 

M 



We show that 



carrying, finite bias situation [1| 
for a wide range of parameters the 2CK fixed point is 
stabilized by a correlation-induced level crossing, forcing 
the degenerate defect levels to be the interacting ground 
states of the defect. In addition to the 2CK Kondo tem- 
perature Tk, the 3LS has another correlation- induced 
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low-energy scale TJ, associated with ground-to-excited- 
state transitions. can be exponentially smaller than 
Tk- We show by explicit renormalization group (RG) cal- 
culations of the non-equilibrium differential conductance 



19| that this can account for conductance spikes at 
elevated bias, much narrower than the ZBA. 

The model. — For the defect levels we employ the 
Abrikosov pseudospin representation, where is the 
creation operator for the defect in state m = 0, ±1, m = 
labeling the (non-interacting) ground state and m — ±1 
the excited doublet [cf. Fig. [T^)] . The defect dynamics 
are subject to the constraint Q — J2m=o±ifmfm — 1- 
The defect coupled to the conduction electron continuum 
is then described by the Hamiltonian, 
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where the first term represents the conduction electron 
kinetic energy, the second one the degenerate local dou- 
blet with level splitting Aq above the defect ground state, 
and the third and fourth terms transitions between the 
local levels induced by conduction electron scattering, 
see Fig. [T^). The conduction electron operators, cj^g-m' 
carry the conserved magnetic spin a — ±1/2 as well as 
an SU(3) index m — 0,±1 which is altered by the defect 
scattering, and which is thought to be associated with 
an electronic orbital degree of freedom centered around 
the defect. The prime in the kinetic energy indicates a 
restricted momentum sum such that X^km = X^k covers 
the complete momentum space. Note that in the rota- 
tional defect model the effective band cut-off reduction 
of Ref. due to screening by high-energy electrons docs 
not occur, because the transitions ni = 1 ^ m — —\ 
are not associated with a charge redistribution. For later 
use we assume that the defect is located in a nanoscopic 
point contact between two metallic leads held at chem- 
ical potentials ii± — ^ ±V/2 [11]. The electron dis- 
tribution function at the location of the defect is then 
= + V/2) + - V/2)]I2, with /o(c^) the 

Fermi distribution function. Throughout, Greek super- 
scripts a, = ± indicate whether an electron is from the 
high or the low potential reservoir. The defect 
SU(3) operators are defined as Sz = flfi — and 



Smn = fmfm ^"^^ the SU(3) operators acting on the elec- 
tronic Fock space, Sz, Smn, are obtained by substituting 

f m ^ X^ 
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in the above expressions. In Eq. ([T]) we 
have chosen a representation of SU(3) which explicitly ex- 
hibits the unbroken symmetry of the SU(2) subgroup in 
the states m = ± 1 . In the following we use dimensionless 
couphngs denoted by the identification N{0)gmn — > ffmn. 

Nonequilihrium RG. — The defect is described by the 
(retarded) impurity propagator Gm(v) = [v — |m|A -|- 
iFm]"^, which, by symmetry, remains diagonal in m, 
even when coupled to the electron system. Here A is 
the renormalized level-spacing (see below), and Tm is the 
interaction-induced decay rate of level m. When a bias 
voltage V is applied, all electrons within at least the volt- 
age window contribute to the scattering. The corrections 
to the electron-defect vertex depend on the energy of the 
incoming electron, oj, and of the in- and out-going defect 
states n, m. To a very good approximation [l8[, the de- 
fect dynamics can be taken on-shell [v = |m|A), so that 
the electron-defect vertex depends on n, m, and on w as 
the only continuous variable. This has led in Ref. [l3| to 
the concept of a perturbative RG for coupling functions, 
parameterized by the energy of the incoming electron, 
i.e. the Hamiltonian, considered at each electron energy 
uj, has its own RG flow. The one-loop RG equations for 
our SU(3) model can be derived as [l8|. 
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FIG. 2: (Color online) Coupling function gl^^^ (ui) for different 
bias voltages V (parameter values as in Fig. [3)|. The peak 
structure of (?qj' (lo) is centered around the renormalized level 
spacing A{D = 0) and is split by the respective bias V. The 
inset shows the defect selfenergy Em(ii'), m — 0,1, calculated 
in bare 2nd order PT. The level spacing renormalization is 
the difference between the two solid dots, SA = Re(Ei(A) — 
So(A)). 
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with g[^^_i = 
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= J± and gmln — for 
j = ±l,m = ±l. The 9 step functfon in Eq. ([2]) controls 
that the RG flow of a particular term is cut off when the 
band cutoff D, symmetrical about the Fermi edge 
in lead 7 = ±, is reduced below the energy of the in- 
termediate scattering electron in that lead. This energy 
is Vlni = w + (|n| - |^|)A(i:)). The Kronecker-(5 factor 
excludes non-logarithmic terms which do not alter the 
impurity state. The energy arguments of the (7's on the 
right-hand-side (RHS) arise from energy conservation at 
each vertex. The exchange terms, exch., on the RHS of 
Eq. (21) are obtained from the direct ones by interchang- 
ing in the g's in- and out-going pseudofermion indices 
and by interchanging m ^ n everywhere. The RHS also 
contains the decay rates [l8| and the renormalized 
level spacing A. Although Ym and A contain no leading 
logarithmic terms, they acquire an RG flow, since they 
are calculated in 2nd order Keldysh perturbation the- 
ory (PT) using the running couplings gmni'jj) 
Specifically, A is obtained during the RG flow as. 
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A = Ao + Rel]i(i/ = A) - YlcY.o{v = 0) , 



(3) 



We find no significant F-dependence of A for 
ln(V^/_Do) 1- Note that the renormalized t/'s have, in 
general, different lo dependence, even if their bare values 
are chosen equal. The typical behavior of the coupling 
functions as solutions of Eq. ([2]) for _D — > is shown in 
Fig. [2] for ,gQ°^ (w) as an example. As expected for reso- 
nant transitions from the excited to the impurity ground 
state, ^0*^^ {uj) has a peak structure centered around the 
renormalized transition energy A. The peak is cut off by 
Tm and split by the applied bias, reflecting the difference 
between the Fermi levels in the two leads. The Kondo 
temperature Tk has been determined here and through- 
out as the cutoff value D for which the coupling reaches 
J_i(w = ±V/2) = 1. 

Equilibrium: 2CK fixed point. — To investigate the 
fixed point structure of the model ([1]), we first analyze 
the RG fiow in equilibrium at temperature T = 0. The 
Greek superscripts have then no relevance. Note that 
even for V = Q the renormalized couplings have an lo- 
dependence. However, in this case it is sufficient to con- 
sider electrons at the Fermi level only, lo — and Vm — 0. 
Due to time reversal symmetry we have in equilibrium, 
gmn{'-^) = gnm{^ + (|"-| " Wc havc chcckcd that 

for V = the results of the RG equations ([2]) coincide 
with those of the familiar equilibrium perturbative RG. 
As the most important point of this paper we find that 
the level spacing A is always renormalized downward by 
the RG flow and, for a wide range of bare parameters 
Jz, Jj_, Ao, and g, crosses the impurity ground level. 
Such a crossing is not forbidden by symmetry, since even 
for completely symmetric initial couplings they have al- 
ready been renormalized to different values when the 
level crossing occurs. Typical RG flow curves A(ln_D) 



are shown in the inset of Fig. [3] It is seen that the level 
crossing occurs generically in the regime ln(D/T/f ) ^ 1, 
where the perturbative RG treatment is still controlled. 
Hence, this behavior can qualitatively be understood by 
2nd order PT as follows (cf. Fig. [21 inset). The imagi- 
nary part of the defect selfenergy Em(t^) is well known 
to have threshold behavior with Iml]m(:/ < 0) = [2l[, 
Because of resonant transitions between the degenerate 
states m = ±1 one has always |ImI]i(^)| > |ImI]o(i^)|, so 
that via Kramers-Kronig the level spacing renormaliza- 
tion in Eq. Re[I]i(A) — 5]o(0)] is necessarily negative 
(the difference between the black dots in Fig. [21 inset). 
We find that even in bare 2nd order PT for Sm(i^) the 
renormalization may be stronger than the bare spacing 
Ao, thus inducing a level crossing. Once the level cross- 
ing has occurred, the system is bound to flow to a stable 
2CK fixed point, since the renormalized ground state of 
the local defect is now doubly degenerate (Kondo degree 
of freedom) , stabilized by spatial parity of the underlying 
crystal lattice, and the channel symmetry of the magnetic 
electron spin is preserved by Kramers degeneracy. Our 
results are summarized in the equilibrium phase diagram 
shown in Fig. [31 The 2CK phase is realized for a sizable 
range of parameters. 

Nonequilibrium differential conductance. — Finally we 
compute the current through an ultrasmall metallic point 
contact in the presence of 2CK defects. We assume that 
a conductance channel is blocked when a 2CK scattering 
from that channel occurs. The current correction due to 
a single 2CK defect is then SI{V) = -f ([(p+ - p_), i/]). 
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We computed this expres- 
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FIG. 3: (Color online) Phase diagram in the — Ao plane 
for Jz/2 — J± = 0.018, bias V = 0, and various values of 
5om ~ flmoi n-jTri = 0, ±1. The symbols mark the line of 
instable fixed points separating the 2CK (left of the line) from 
the potential scattering (right of the line) low-temperature 
phase. The solid lines are straight line fits to the vs. Ao 
plots and show that the level crossing is essentially an 0(j\) 
PT eff'ect. The inset shows typical RG flow diagrams of the 
level spacing A(D) in a) the 2CK phase (level crossing) and 
b) the potential scattering phase (no level crossing). 
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the model suggests a unified microscopic explanation for 
the essential features of the conductance anomalies ob- 
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FIG. 4: DifFerential conductance correction dl/dV per defect 
for bias V near the renormalized level spacing A{D = 0), 

> « IQ-^Do. J. 12 = Jx = 0.016, g^l' = g^^ = 0.6Jx, 
Ao = STk, A(D = 0) = -1.82 Tk. Right inset: d//dV 
for a larger bia s ra nge. The dashed line is a sketch of the 
2CK- induced y|F|-ZBA of width Tk, with the numerically 
determined Tk, as indicated. Left inset: Experimental d//dl/ 
curve taken from Ref. fTl| for comparison. 



sion for T = 0, as outlined in Ref. by evaluating 
the 2nd order Keldysh diagram shown in Fig. [T}d) [l^ , 
using the renormalized coupling functions obtained from 
Eq. ([2]), and obtained the voltage dependent nonequi- 
librium conductance dl /dV . Details of the lengthy cal- 
culation will be presented elsewhere. The results are 
displayed in Fig. IH showing conductance spikes much 
narrower than Tk at a bias V corresponding to the 
renormalized level spacing for that bias, V = ±|A|. 
These spikes, with a steep step on the low-V side and 
a wider (logarithmic) decay on the high-T^ side, are anal- 
ogous to the conductanc e p eaks of a Kondo quantum 
dot in a magnetic field [l8|. They are due to transi- 
tions between the (renormalized) ground- and excited- 
state defect levels, enhanced by Kondo scattering. De- 
spite their complicated shape arising from a combina- 
tion of many terms, their overall narrow width can be 
understood from the fact that the model ^ contains 
at least one more low-energy scale in addition to Tk'- 
While Tk is to leading log order determined by fluctu- 
ations between the degenerate defect levels (m = ±1), 
i.e. by Jj_, Tk — -Doexp[— l/4Jj_], the conductance has 
in 3rd order [cf. Fig. [T}:)] a logarithmic divergence at 



V 



±|A|, which is governed by the couplings g^Q^ 
dl/dV ^ 45^;) ln[e(|F| - |A|)(|y| - m/D,]. 
The width of this peak, although additionally cut off by 
the nonequilibrium rates F™, may thus be estimated as 



T* 



i.e. for f/j^Q < Jl it is exponcn- 
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tially smaller than Tk 

In conclusion, we have shown that a three level defect 
with partially broken SU(3) symmetry and degenerate 
excited states in a metal has a sizable, stable 2CK fixed 
point phase and that this model can describe conduc- 
tance spikes exponentially more narrow than Tk- Hence, 



served by Ralph and Buhrman [11| . The reduction of the 
ZBA and the splitting of the spikes in a magnetic field 
B [lH can be explained by a coupling of B to the defect 
angular momentum and the resulting Zeeman splitting 
of the levels m = ±1. This will be discussed elsewhere. 
The size of the anomalies ~ 100 /h might be due to 
a corresponding number of 2CK defects in the contact. 
The very sharp experimental distribution of Tk UJ may 
occur, because the rotational defects form at equivalent 
interstitial sites in a (nearly perfect) Cu lattice or because 
only those defects with sufficiently high Tk are observ- 
able [13]. — We wish to thank J. von Delft, A. Rosch, 
and especially A. Zawadowski for numerous useful discus- 
sions. This work was supported by DFG through SFB 
608 and grant No. KR1726/1. 
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